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Abstract 

We detected broadband coherent terahertz (THz) emission from multiferroic BiFe03 after illu- 
minating a high-quality bulk single ferroelectric domain crystal with a ~100 fs optical pulse. The 
dependence of the emitted THz waveform on the energy and polarization of the optical pulse is 
consistent with the optical rectification mechanism of THz emission. The THz emission provides 
a sensitive probe of the electric polarization state of BiFeOs, enabling applications in ferroelec- 
tric memories and ferroelectric domain imaging. We also report room-temperature THz optical 
constants of BiFeOs • 
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Multiferroic BiFeC>3 (BFO) offers the possibility of electrical manipulation of its mag- 
netic state, which can be used in devices such as electrically-controlled magnetic memory 
elements^. At room temperature, single- crystalline BFO is found in a G-type antiferromag- 
netic state (Xjv=640 K) with a modulated cycloidal spin structure with a long modulation 
period of ~620 A-^£. The material is also a robust ferroelectric with a high Curie tem- 
perature Tc of 1,103 K and large ferroelectric polarization of ~100 /iC/cm 2 -. Nonlinear 
optics plays an important role in characterization of ferroelectrics, as the lack of inversion 
symmetry in ferroelectric crystals results in the presence of strong second-order nonlinear- 
ities; nonlinear optics is also an important area for ferroelectric applications, along with 
transducers and data storage. Linear and nonlinear optical properties of BFO have been 
reported^&i&ii. Recently, it has been proposed that terahertz (THz) radiation from BFO 
thin films can be used in ferroelectric domain imaging and optical readout of the state of fer- 
roelectric memories^ 1 ^ 1 ^. In that work, the THz wave was generated by ulrafast modulation 
of ferroelectric polarization when charge carriers were injected in the film by a femtosecond 
laser pulse. 

In this Letter, we report coherent generation of THz waveforms by femtosecond 800-nm 
optical pulses from ferroelectric BFO in reflection and transmission geometries and also 
provide the THz- frequency optical constants of BFO. The THz wave is generated via optical 
rectification that relies on the second-order nonlinearity. The amplitude of the THz wave 
detected in the time-domain depends on the direction of the ferroelectric dipole in the BFO 
crystal. In contrast, a measurement of the intensity of the second-harmonic light does 
not provide this information^ 1 ^. Thus, THz emission provides a sensitive probe of the 
polarization state of the BFO crystals. The presented work was performed on a high-quality 
single ferroelectric domain crystal of BFO, which is evidenced by the pronounced polarization 
dependence of both the optical rectification and the measured THz optical constants. Such 
high-quality crystals are indispensable for the reliable characterization of the basic material 
properties' 1 ^. 

For this work, bulk single ferroelectric domain BFO crystal was grown using a Bi203 flux. 
The ferroelectric phase possesses the rhombohedral R3c structure that can be obtained 
by a slight elongation of the pseudocubic perovskite unit cell along the (111) direction. 
Our measurements were performed at room temperature on a 0.22 mm thick (001) cu bj C - 
oriented crystal mounted on a 3-mm diameter aperture. Coherent THz transients were 
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generated upon illumination of the crystal by 800-nm femtosecond pump pulses provided 
by Ti:Sapphire regenerative amplifier operating at 1 kHz and focused to a 2-mm diameter 
spot. The maximum pump power of 150 mW used in these measurements corresponds to 
the fluence of 4.8 mJ/cm 2 . The THz transients emitted from BFO were focused on a [110]- 
oriented ZnTe crystal for coherent detection via electro-optic sampling^^. The THz optical 
constants of the BFO crystal were measured in transmission in a time-domain spectroscopy 
system that uses photoconductive antennas for THz generation and detection 19 . 

Figure [1] shows the geometry of THz generation measurement in BFO. The incident 
femtosecond pump pulses are polarized in the horizontal x-z plane. The ferroelectric dipole 
moment in BFO points along the long diagonal of the distorted perovskite structure and 
makes a 35° angle with the surface of our crystal. Figure [T] displays the THz transients 
emitted from the crystal in transmission and reflection geometries. Both in transmission 
and reflection, the main THz pulse is followed by the Fabry-Perot (or etalon) reflection 
within the sample. 

When a lightwave of frequency to is incident on a crystal with a second-order optical 
nonlinearity the polarization induced by the lightwave contains the nonlinear part 
p nl = x (2)( cosa;t )2 = i/2 x ( 2 )(l + cos2cj) = P NL (0) + P NL (2u). The latter term in 
the last expression describes the second harmonic generation. The former term describes 
a zero-frequency polarization, which is truly a constant (DC) polarization induced in the 
crystal when the incident lightwave is a sinusoidal field of infinite duration. Induction of DC 
polarization due to the second-order nonlinearity is known as optical rectification. When 
the incident lightwave is a 100-fs-long pulse (Fig. UK), the induced P NL (0) only exists for the 
duration of the pulse and represents a time- varying electric dipole P(t). This time-dependent 
dipole emits a THz electromagnetic wave E T Hz(t) ~ d 2 P(t)/dt 2 . To confirm that optical 
rectification is responsible for the observed THz emission, we measure the dependence of the 
emitted THz amplitude on the angle between the [110] cu 6j c crystallographic direction and 
the y axis while rotating the crystal about its normal (Fig. UK). The angular dependence 
of the emitted THz field is consistent with the optical rectification mechanism. We model 
the induced second order THz-frequency polarization Pj = xfjkEj(uj)E k (uj) using the set 
of nonlinear optical constants of BFO films found by Kumar et a/ 10 . The emitted THz 
field E THz {t) has the same angular dependence as the induced P{t). The solid line in 
Fig. [2] represents the calculated induced polarization and shows excellent agreement with 
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the measured THz field amplitude. The meaning of the negative THz amplitude is illustrated 
in the inset of Fig [2j The optical rectification mechanism of the THz emission is supported 
further by the linear dependence of the THz field amplitude on the pump power (inset of 

Fig. ED- 

The ferroelectric axis of BFO can be along any of the four long diagonals of the pseudocu- 
bic unit cell (Fig. Hk). Correspondingly, eight different ferroelectric domains are possible in 
the crystal. The pronounced angular dependence of THz emission (Fig. [2]) demonstrates that 
our crystal predominantly consists of a single ferroelectric domain, which is a result of crys- 
tal growth below the Curie temperature. The single-domain nature of these millimeter-sized 
crystals was confirmed independently by neutron diffraction^. 

A different mechanism for THz emission from BFO thin films was described by Takahashi 
et alM^^. In their work, the emission is due to the injection of photoexcited charge carriers 
and ultrafast modulation of spontaneous polarization Pf. In our measurements, similar 
processes involving the injection of photocarriers can be ruled out because the energy of our 
800 nm, 1.5 eV pump pulses is lower than the direct optical gap of BFO. To confirm this, 
we carried out optical-pump THz-probe study of the BFO crystal, in which no changes to 
the THz optical constants were induced by the 800-nm pump pulses. The absence of THz 
conductivity modulation indicates that no free carriers are created by the 800-nm lights. In 
contrast, Takahashi et al used 400-nm pump light in their work to inject the photocarriers. 
Our BFO crystal is opaque to the 800-nm light, even though this photon energy is lower 
than the measured bandgap. The residual absorption is most likely due to a certain level of 
present defects. The same defects are expected to act as traps for the few carries that are 
created by the absorption. Thus, the residual absorption in our BFO crystal does not result 
in a measurable change in THz conductivity. 

We determine the spectral content of THz pulses emitted from BFO by taking a Fourier 
transform of the measured time-domain transients and compare it to THz emission from 
0.5-mm thick [110]-oriented ZnTe crystal (Fig. [3]). The THz spectrum emitted from BFO in 
reflection is very similar to that of the emission from ZnTe, although the amplitude of ZnTe 
emission is approximately 25 times larger than that of BFO at the same pump power. The 
differing emission conditions in ZnTe and BFO, i.e., different crystal thicknesses, emission 
direction and phase matching, suggest that in both cases the ~ 2.5 THz bandwidth may 
be limited by the frequency response of the electro-optic-sampling THz detector used in the 
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measurement. 

The spectrum of the THz pulse emitted from BFO in transmission is considerably nar- 
rower than the spectrum of the pulse emitted in reflection (Fig. [3J in this Figure all spectra 
are normalized by the maximum amplitude of the ZnTe spectrum). We attribute this to 
the absorption of the higher frequencies in the THz pulse by BFO as the pulse propagates 
through the crystal. We point out that even though the optical absorption at 800 nm is 
negligible compared to the above-bandgap absorption in BFO, the residual absorption is 
sufficiently strong to prevent the 800-nm light from propagating far into the BFO crystal. 
We estimate that the penetration depth of the 800-nm light is about 10 /im, which is only a 
small fraction of the total 220 /im thickness of our BFO crystal. (Our BFO crystal is com- 
pletely opaque to the eye. A different BFO crystal became transparent after we polished 
it down to 10 /j,m thickness.) Therefore, the THz wave is emitted in a thin surface layer of 
our BFO crystal and then propagates through and is absorbed by the remaining portion of 
the crystal. To confirm this, we measured THz absorption in BFO (upper panel of Fig. H]) 
and found that the crystal is indeed opaque above ~ 1.6 THz for the light polarized along 
the ferroelectric polarization Pf. The THz optical constants extracted from the measure- 
ment (lower panel of Fig. H|) allow an estimate of the THz pulse amplitude emitted from the 
10- /xm-thick layer in which the pump pulse is absorbed. Thus calculated THz amplitude 
is comparable to the measured amplitude of the THz pulse emitted in reflection. The two 
amplitudes are expected to be the same because phase matching is unimportant in a suffi- 
ciently thin crystal. More evidence for THz emission occurring in a thin surface layer of the 
crystal is provided by the broadened Fabry-Perot reflection of the main THz pulse detected 
in reflection (Fig. [1]). The time-domain width of the Fabry-Perot pulse is the same as that 
of the main transmitted THz pulse, which means that the main reflected pulse is narrow 
in the time domain and has a much broader spectrum in the frequency domain because it 
does not traverse an appreciable thickness of BFO material and therefore does not lose the 
high-frequency components in its spectrum. 

In conclusion, we demonstrate the emission of coherent THz radiation from BFO due 
to optical rectification of ~ 100 fs 800-nm pump pulses both in reflection and transmis- 
sion geometries. This mechanism of THz emission is different from the previously reported 
emission via ultrafast modulation of the spontaneous ferroelectric polarization^^^. The 
potential applications of coherent THz emission include ferroelectric domain imaging and 
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non-destructive readout of the state of ferroelectric memory elements. Since optical recti- 
fication does not require the absorption of the pump pulse, our demonstration opens the 
door for future studies using pump wavelengths longer that the 442-nm bandgap of BFO, in- 
cluding the important 1.3 /im and 1.55 /xm telecommunications wavelengths. We also report 
room-temperature THz- frequency optical constants of BFO, which will be indispensable in 
modeling THz emission and propagation in those devices. 
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FIG. 1: (Color online) (a) Geometry of THz generation in (001) cu fej c -oriented BFO. The pump light 
is polarized in the x-z plane. Static ferroelectric polarization P is along the [111] direction (the 
diagonal of the cube), (b) THz waveforms detected in transmission and reflection geometry, pump 
fluence 4.8 mJ/cm 2 . The main THz pulse and the first Fabry-Perot reflection are shown. 
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FIG. 2: (Color online) Amplitude of the main THz pulse as a function of the angle cj> between 
the y axis and the [110] direction in the crystal while rotating the crystal about the z axis. The 
solid line shows the simulated nonlinear polarization P(t). The inset shows the main THz pulse at 
4> = 270° and <fi = 90°, the positions of maximum positive and negative THz amplitudes. 
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FIG. 3: (Color online) Frequency-domain spectra of the waveforms emitted from BiFeOs and from 
the [110] oriented ZnTe crystal. The THz waveforms were measured with the pump power of 150 
mW and 28 mW, respectively. The inset shows the amplitude of the THz wave from BiFeO"3 
detected in transmission as a function of pump power. 
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FIG. 4: (Color online) Upper panel: Amplitude transmission spectra of BiFeO"3 crystal for THz 
light polarized in the plane of the electric dipole (solid line) and perpendicular to the plane (dotted 
line) . Lower panel: real and imaginary parts of the refractive index for the polarizations displayed 
in the upper panel. 
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